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1. Introduction 

In this paper we are concerned with the coefficients in the three-term recurrence 
relations for semi-classical orthonormal polynomials, specifically for generalizations of 
the Charlier and Meixner polynomials which are discrete orthonormal polynomials. It 
is shown that these recurrence coefficients for the generalized Charlier polynomials and 
generalized Meixner polynomials can respectively be expressed in terms of Wronskians 
that arise in the description of special function solutions of the third Painleve equation 
(Pm) 

d 2 w 1 fdw\ 2 ldw Aw 2 + B „ , D , . 

tt = - hr -~A- + + c w 3 + -, i.i 

dz w \dz ) z dz z w 

where A, B, C and D are arbitrary constants, and the fifth Painleve equation (P v ) 
d 2 w _ ( 1 + 1 \/dwA 2 _ ldw + (w- l ? ( Aw+ B 



dz \2w w — 1 / \ dz J z dz z 2 \ w 

+ Cw + P W ( W + l) 
z w — 1 

Wronskians for special function solutions of Pni are expressed in terms of modified 
Bessel functions and for Py in terms of confluent hypergeometric functions. 

The relationship between semi-classical orthogonal polynomials and integrable 
equations dates back to the work of Shohat [55] in 1939 and later Freud [29] in 1976. 
However it was not until the work of Fokas, Its and Kitaev [23, 24] in the early 1990s 
that these equations were identified as discrete Painleve equations. The relationship 
between semi-classical orthogonal polynomials and the (continuous) Painleve equations 
was demonstrated by Magnus [40] in 1995. A motivation for this work is that recently 
it has been shown that recurrence coefficients for several semi-classical orthogonal 
polynomials can be expressed in terms of solutions of Painleve equations, see, for 
example, [2, 3, 5, 6, 8, 10, 11, 12, 13, ?, 16, 17, 19, 20, 21, 26, 28, 56, 60, 61]. 

This paper is organized as follows: in §2 we review properties of Pm (1.1) and 
Pv (1-2), including special function solutions and the Hamiltonian structure of these 
equations; in §3 we review properties of orthogonal polynomials and discrete orthogonal 
polynomials; in §4 we derive expressions for the recurrence coefficients for the generalized 
Charlier polynomials in terms of Wronskians that arise in the description of special 
function solutions of Pm; in §5 we derive expressions for the recurrence coefficients for 
the generalized Meixner polynomials in terms of Wronskians that arise in the description 
of special function solutions of P v ; and in §6 we discuss our results. 

2. Painleve equations 



The six Painleve equations (Pi-Pvi) were first discovered by Painleve, Gambier and 
their colleagues in an investigation of which second order ordinary differential equations 
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of the form 

d 2 w ^ / dw \ .„ _ . 

where F is rational in dw/dz and u> and analytic in z, have the property that their 
solutions have no movable branch points. They showed that there were fifty canonical 
equations of the form (2.1) with this property, now known as the Painleve property. 
Further Painleve, Gambier and their colleagues showed that of these fifty equations, 
forty-four can be reduced to linear equations, solved in terms of elliptic functions, or 
are reducible to one of six new nonlinear ordinary differential equations that define new 
transcendental functions (see Ince [34]). The Painleve equations can be thought of as 
nonlinear analogues of the classical special functions [15, 22, 31, 36], and arise in a wide 
variety of applications, for example random matrices, see [25, 54] and the references 
therein. 

The Painleve equations Pn-Pyi possess hierarchies of solutions expressible in terms 
of classical special functions, cf. [15, 31, 43] and the references therein. For P m (1.1) 
these are expressed in terms of Bessel functions [45, 46, 52], which we discuss in §2.1, 
and for P v (1.2) in terms of confluent hypergeometric functions (equivalently, Kummer 
functions or Whittaker functions) [51, 43, 62], which we discuss in §2.3. 

Each of the Painleve equations Pj-Pvi can be written as a (non-autonomous) 
Hamiltonian system. For P n i (1.1) and P v (1.2) these have the form 

dg dUj dp dH 3 

Z T~ = ' Z T~ = — a i J = m > v 
dz op dz oq 

for a suitable Hamiltonian function T-Lj(q,p,z) [37, 49, 51, 52], which we discuss in 
§2.2 and §2.4. Further, the function o~(z) = T-Lj(q,p,z) satisfies a second-order, second- 
degree equation, which is often called the Jimbo-Miwa- Okamoto equation or Painleve a - 
equation, whose solution is expressible in terms of the solution of the associated Painleve 
equation [37, 50, 51, 52]. Hence there are special function solutions of these equations, 
which are also discussed in §2.2 and §2.4. 

2.1. Special functions solutions of the third Painleve equation. 

In the generic case when CD ^ in P m (1-1), then we set C = 1 and D = —1, without 
loss of generality, so in the sequel we consider the equation 

d 2 w 1 (dw\ 2 Idw Aw 2 + B , 1 , n n . 

— = -( — ) --— + + w 3 --. (2.2) 

dz w \ dz / z dz z w 

Special function solutions of (2.2) are expressed in terms of Bessel functions, see 
[45, 46, 52]. 

Theorem 2.1. Equation (2.2) has solutions expressible in terms of Bessel functions if 
and only if 

eiA + e 2 B = 4n + 2, (2.3) 
with n G Z and S\ — ±1, £2 = ±1 independently. 
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Proof. See Gromak [30], Mansfield and Webster [41] and Umemura and Watanabe [59]; 
also [31, §35]. □ 

An alternative form of P m , due to Okamoto [49, 50, 52], is obtained by making the 
transformation w(z) = u(t)/y/i, with t = ^z 2 , in (2.2) giving 

fu_l/du\ 2 ldu Au^_ B v? 1 

d? ~u\dt) ~lltt + l>¥ + 2l + l 2 ~u' 
which is known as Pm'. Equation (2.4) has solutions expressible in terms of solutions 
of the Riccati equation 

t— = e x u + vu + e 2 t, (2.5) 

if and only if A and £> satisfy (2.3). To solve (2.5), we make the transformation 

d 

u(t) = -Eit—lnip v (t), 
then ip v (t) satisfies 

which has solution 



t-nr + (1 - v)-TT + eie*^ = ( 2 - 6 ) 



' f' 2 {dJ v (2Vi) + C 2 Y v (2Vi)} , if ei = l, e 2 = l, 

r"/ 2 {d J^V*) + (72^(2^)}, if ei = -l, e 2 = -l, 

t" /2 {^i/,(2v / t) + C , 2^(2 v / t)}, if ei = l, £ 2 = -l, 

^r^ 2 {c 1 / ! ,(2v / t) + C2^(2v / t)}, if £i = -l, £ 2 = 1, 



(2.7) 



with C\ and C 2 arbitrary constants, and where J u (z), Y u (z), I v {z) and K v [z) are Bessel 
functions. 

2.2. Hamiltonian structure for the third Painleve equation. 
The Hamiltonian associated with Pnr (2.4) is 

n m >(q,P,t) = qV- {q 2 + Aq - t) p + ±{A + $oo)q, (2.8) 
with i? an d "^oo parameters, where p and g satisfy 



^ = 2g 2 p-g 2 -^ g + t, 



t^j = -2qp 2 + 2qp + $ p - |(^ + #oo), 



(2.9) 



L dp 

'dt 

see Okamoto [49, 50, 52]. Eliminating p in (2.9) then q = u satisfies Pnr (2.4) with 
parameters (A, B) = (— 2-# 00 , 2(-# + !))• Eliminating q in (2.9) then p satisfies 



d 2 p 1^1 1 \ f dp^ 1 dp 2p(p — 1) 



dt 2 2 \p p - 1 J \ dt J t dt t 

.'W.fcM.itMj, ( , 10) 

8t 2 { p p — 1 1 
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Making the transformation p(t) — 1/[1 — w(z)\, with z — t, in (2.10) yields 

/ d-uA 2 



d 2 w 
d? 



1 1 

+ 



— 



1 dw 

z dz 



2w w — 1 / \ dz J 

<W-l) 2 /(^o+^oo) 2 ^ (^0-^oo) 2 



+ 



2w 

; 

Z 



8u, ■ -■ < 2 ' U > 

which is P v (1.2) with parameters (A, B, C, D) = (§(tf + ^oo) 2 , -|(#o - $oo) 2 , -2, 0), 
illustrating the well-known relationship between Pni and Py with 5 = 0, cf. [31, §34]. 

The second-order, second-degree equation satisfied by the Hamiltonian function is 
given in the following theorem. 

Theorem 2.2. The Hamiltonian function 

a(t; O , = tUw, (q, p, t)-\t+ 

with Uni' (q,p,t) given by (2.8), satisfies the second-order, second-degree equation 



dV 
dt 2 



+ <4 



(da 



\dt 



- 1 



da 
dt 



da 
dt 



(2.12) 



Conversely, if a(t) satisfies (2.12) then the solution of the Hamiltonian system (2.9) is 
given by 



d 2 a n n da 
t— j - 2$ — + $ a 
,_dr dt 



da , 

pff) = w + i 



Proof. See Okamoto [50, 52]; see also Forrester and Witte [27]. 

The special function solutions of (2.12) are given in the following theorem. 
Theorem 2.3. Let r njI/ (t) be the determinant given by 

w.) ••• sr 1 ^) 

SM 5 2 {^) ••• <w„) 



(2.13) 



□ 



St. = t 



dr 



(2.14) 



5r\^) w„) ■■■ 5 2n - 2 (w) 

with ipuit) a solution of (2.6). Then special function solutions of (2.12) are given by 



a(t) = t— (In r ntV (t)) + \e x e 2 t + \v 2 + \n{\ - e x v) - \n 2 , 
/or t/ie parameters (i?o,$oo) = (i^ + n, e^C^ — 



(2.15) 



Proof. See Okamoto [52]; see also Forrester and Witte [27]. □ 
The determinant r„^(t) given by (2.14) is often called a "r-function" , see [52]. 
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2. 3. Special functions solutions of the fifth Painleve equation. 

In the generic case when D ^ in Py (1-2), then we set D = — |, without loss of 
generality, so in the sequel we consider the equation 

d 2 w _ ( 1 1 \/dwA 2 \dw + (w- l) 2 / B 



dz \2w w — 1 / \ dz J z dz z 2 \ w 

Cw_ W ( W + l) 
z 2(w-l) K J 

Special function solutions of (2.16) are expressed in terms of confluent hypergeometric 
functions (equivalently Kummer functions or Whittaker functions), see [51, 43, 62]. 

Theorem 2.4. Equation (2.16) then has solutions expressible in terms of Kummer 
functions if and only if 

a + b + e 3 C = 2N + 1, (2.17) 

or 

(a - N)(b - N) = 0, (2.18) 
where N G N, a = E\\f2A and b = 62V— 2B, with Ej = ±1, j = 1, 2, 3, independently. 

Proof. See Okamoto [51], Masuda [43] and Watanabe [62]; also [31, §40]. □ 

In the case when iV = in (2.17), then (2.16) has solutions in terms of the associated 
Riccati equation 

z^- = aw 2 + (b- a + e 3 z)w -b. (2.19) 
If a 7^ 0, then (2.19) has the solution 

w(z) = ---^-]n<p(z), 
a dz 



where ip(z) satisfies 

d 2 dip 

z 2 — % + z(a — b + 1 — £32)- ab(p = 0, 

dz dz 

which has solution 

z b {dM(b, l + a + b,z) + C 2 U(b, 1 + a + b, z)} , if e 3 = 1, 

z b e~ z {dM(l + a, 1 + a + 6, z) + C 2 U(1 + a, 1 + a + b, z)} , if e 3 = -1, 

with C\ and C 2 arbitrary constants, and where M(a,j3,z) and U(a,j3,z) are Kummer 
functions. 



(fi(z) 
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2.4- Hamiltonian structure for the fifth Painleve equation. 
The Hamiltonian associated with (2.16) is 

Hy(q,p, z) = q(q - 1) V - {(6 + d)q 2 - (26 + 7? - z)q + b}p - \{a 2 - (b + d) 2 }q, (2.20) 
with a, b and d parameters, where p and q satisfy 

zp- = 2q(q - l) 2 p - (6 + $)q 2 + (26 + $ - z)q - 6, 

z-f = -(3q - l)(g - l)p 2 + 2(6 + $)qp - (26 + - z)p + ±{a 2 - (6 + tf) 2 }, 

see Jimbo and Miwa [37] and Okamoto [49, 50, 51]. Eliminating p then q = w satisfies 
(2.16) with parameters (A,B,C) = (±a 2 , -|6 2 , -d - 1). 

The second-order, second-degree equation satisfied by the Hamiltonian function is 
given in the following theorem. 

Theorem 2.5. The Hamiltonian function 

a(z; a, 6, 0) = ft v (g,P, *) + K 26 + d ) z ~ K 2b + ^> ( 2 - 22 ) 
with T-Ly(q,p, z) given by (2.20), satisfies the second- order, second-degree equation 



d ; 



(7 



d^ 



-<2 



dcr\ do" 



da 

d7 



+ Kj)= 0, (2.23) 



/or £6e parameters 

K 1 = \d + \a, K2 = \ti-\a, K Z = -\d-\b, n A = -\<8 + \b. (2.24) 

Conversely, if a{z\ a, 6,1?) satisfies (2.23) then the solution of the Hamiltonian system 
(2.21) is given by 



q(z) = 



za 



2(a') 2 -za' + a 



p{z) = 



II n/ l\2 i / 

2(7 — ) -|- zo~ — a 



2{a' + \$ -\a){a' 
where ' = d/dz. 

Proof. See Jimbo and Miwa [37] and Okamoto [50, 51]. 

The special function solutions of (2.23) are given in the following theorem. 
Theorem 2.6. Let W n (ip) be the determinant given by 

5 2 M ... saw 



(2.25) 



□ 



and define (p a ,p(z) by 

<p a A z ) = CiM(a, (3, z) + C 2 U(a, p, z), 



A - d 

5 > = z Tz> 



(2.26) 



Discrete orthonormal polynomials and the Painleve equations 8 

with Ci and C 2 arbitrary constants, and where M(a,j3,z) and U(a,(3,z) are Kummer 
functions. Then special function solutions of (2.23) are given by 

<j(z) = z-^{lnW n (ip a ,p)) - J(3n + 2a-(3-l)z 

- §n 2 - \(2a - 3/3 - l)n - \{2a - (3 - l) 2 , (2.27a) 

a{z) = ^(ln W„(A«,/0) " K 3 ™ + 2a - /3 - l)z 

- §n 2 - i(2a + /3 - l)n - |(2a - f3 - l) 2 , (2.276) 

/or the parameters 

Kl = \(2a-0-n-l), k 3 = j(2a-/3 + 3n- 1), 

/t 2 = -i(2a + /3 + n - 3), « 4 = -f (2a - 3/3 + n + 1), 

and 

<r(z) = ^(ln^f/e-^)) + I(3n -2a + f3-l)z 

- |n 2 + i(2a -3/3 + l)n - \{2a - f3 + l) 2 , (2.28a) 

/or the parameters 

K! = \(2a-p + n + l), k 3 = i(2a-/3-3ra + l), 

K2 = - J(2a + /3 - n - 1), k 4 = -J(2a - 3/3 - n + 3). 



(2.27c) 



(2.286) 



Proof. This result can be inferred from the work of Forrester and Witte [27] and 
Okamoto [51]. □ 

3. Orthonormal polynomials 

3.1. Continuous orthonormal polynomials 

Let p n (x), for n G N, be the orthonormal polynomial of degree n in x with respect to a 
positive weight u>(x) on (a, b), which a finite or infinite open interval in R, such that 

Pm(^)Pn(^) ^(x) dx = S min , 



I 

J a 



with 5 m ^ n the Kronekar delta. One of the most important properties of orthogonal 
polynomials is that they satisfy a three-term recurrence relationship of the form 

xp n (x) = a n+1 p n+1 (x) + b n p n (x) + a n p n _i(x), (3.1) 

where the coefficients a n and b n are given by the integrals 

pb rb 

a n = xp n (x)p n ^ 1 (x)u(x)dx, b n = xp\ (x) u(x) dx, 

J a J a 

with p-i(x) = 0. The coefficients in the recurrence relationship (3.1) can also be 
expressed in terms of determinants whose coefficients are given in terms of the moments 
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associated with the weight uj(x). Specifically, the coefficients a n and b n in the recurrence 
relation (3.1) are given by 

2 _ An+iA n _i ^ A n+1 ^ An 



(3.2) 



where the determinants A n and A n are given by 



Mo Mi 

Ml M2 

Mn— 1 Mn 



Mn-1 
M« 

M2n-2 



Mo Ml 
Mi M2 

Mn— 1 M« 



Mrt-2 Mra 
Mn-1 Mn+1 

M2n-3 M2n-1 



(3.3) 



for n = 0, 1, 2, . . . , with A = 1, A_i = and A = 0, and the fcth moment, is given 
by the integral 



rb 

Mfe = / ic fc w(^) dx. 

J a 



(3.4) 



A characterization of classical orthogonal polynomials (such as Hermite, Laguerre 
and Jacobi polynomials), is that their weights satisfy the Pearson equation 



-^[a(x)u(x)] = t(x)u(x), 



(3.5) 



where o~(x) is a monic polynomial with deg(a) < 2 and t(x) is a polynomial with 
deg(r) = 1, cf. [1, 7, 14]. If the weight function u(x) satisfies the Pearson equation 
(3.5) with either deg(a) > 2 or deg(r) > 1, then the orthogonal polynomial is said to 
be semi- classical, cf. [32, 42]. 

For further information about orthogonal polynomials see, for example, the books 
by Chihara [14], Ismail [35] and Szego [58]. 



3.2. Discrete orthonormal polynomials 

One can also define orthogonal polynomials on an equidistant lattice, rather than an 
interval. Consider the discrete orthonormal polynomials {p n ( x )}, n — 0, 1, 2, . . . , with 
respect to a discrete weight oo(k) on the lattice N 

oo 

^Pn{k)p m (k)uj(k) = 5 

k=0 

which also satisfy the recurrence relation (3.1). 

The moments fx n of the discrete weight u(k) are given by 

oo 

fi n = ^k n uj(k), 71 = 0,1,2,... , 

k=0 

and, as for the continuous orthonormal polynomials in §3.1 above, the coefficients in the 
recurrence relation are given by (3.2), with the determinants A n and A n given by (3.3). 
In the special case when the discrete weight has the special form 

u(k) = c(k)z k , z > 0, 
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which is the case for the Charlier polynomials C n (k; z) and the Meixner polynomials 
M n (k;a,z) (see §4.1 and §5.1 below, respectively), then 



^(z) = J2 k n c{k)z k = s:m, 8 z = z—. 



k=0 



(3.6) 



Consequently the determinants A n (z) and A n (z) given by (3.3) have the form 
A n (z) = 



fJ>o 


4(/io) • 








%M ■ 


■ 


J 














■ sr 2 M 






8 2 z (no) ■ 








5 n M . 


■ e- 3 w 





A n (z) 



respectively. Hence we have the following result. 

Theorem 3.1. If the moment /i n (z) has the form (3.6), then the determinants A n (z) 
and A n (z) can be written in the form 

A n (z) = W n M, A n (z) = 5 z W n M, (3.7) 

where W n (ip) is defined by (2.26). 

Some properties of the recurrence coefficients a n (z) and b n (z) are given in the 
following theorem. 

Theorem 3.2. If the moment /i n (z) has the form (3.6), then the recurrence coefficients 
a n (z) and b n (z) in (3.1) satisfy the Toda system 

$z(a 2 n ) = aKbn-bn^), 5 z (b n ) = a 2 n+1 - a 2 n . (3.8) 

Further the determinant A n (z) satisfies the Toda equation 

A n+ iA n _i 



£(lnA n ) 



A 2 



n = l,2,... 



Proof. See [38, 52, 57]; see also [47]. 



□ 



Theorem 3.3. The recurrence coefficients a n (z) and b n (z) in (3.1) can be expressed in 
the form 

W„+i(/i 



a 2 n (z) = 5 2 z (\n b n (z) = <U In 



Proof. Applying Theorem 3.1 to (3.2) gives the result. 



(3.9) 
□ 
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As discussed in §3.1 above, classical orthogonal polynomials are characterized by the 
Pearson equation (3.5). Analogously discrete orthogonal polynomials are characterized 
by the discrete Pearson equation 



a(k) is a monic polynomial with deg(cr) < 2 and r{k) is a polynomial with deg(r) = 1. 
A discrete Pearson equation can also be defined using the backward difference operator 



If the discrete weight u(k) satisfies (3.10) with either deg(a) > 2 or deg(r) > 1, then 
the discrete orthogonal polynomial is said to be semi-classical, cf. [18]. 

For further information about discrete orthogonal polynomials see, for example, the 
books by Beals and Wong [4, Chapter 5], Chihara [14, Chapter VI], Ismail [35, Chapter 
6] and Nikiforov, Suslov and Uvarov [48]. 

4. Charlier polynomials and generalized Charlier polynomials 

4-1. Charlier polynomials 

The Charlier polynomials C n (k; z) are a family of orthogonal polynomials introduced in 
1905 by Charlier [9] given by 



where 2-?o( a ; b; ; z) is the hypergeometric function and L„ (z) is the associated Laguerre 
polynomial, see [4, 14, 35, 53]. The Charlier polynomials are orthogonal on the lattice 
N with respect to the Poisson distribution 



A[a(k)u(k)] = r(k)uj(k), 
where A is the forward difference operator 

A/(fc) = /(*+!) -/(*), 



(3.10) 



V/(fc) = /(*)-/(*-!). 



C n (k; z) = 2 F (-n, -k; ; -l/z) = {-l) n n\L^ (-1/z) , 



z>0, (4.1) 



and satisfy the orthogonality condition 



(4.2) 




k=0 



The weight (4.2) satisfies the discrete Pearson equation (3.10) with 

a(k) = k, r(k) = z — k. 
From (4.2), the moment (J>o(z) is given by 



Hence from Theorem 3.1, the Hankel determinant A n (z) is given by 



n-l 



A n (z) = W n (/io) = «»(»- 1 )/ 2 e" n^ ! ), 



fc=i 
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and so from Theorem 3.3 the recurrence coefficients are given by 

a 2 n (z) = 5 2 z (\nW n (i2 ))=nz, b n (z) = 5 Z ( In ^TT^- ) = n + z. 



W„(/i ) 



4-2. Generalized Charlier polynomials 



Smet and van Assche [56] generalized the Charlier weight (4.2) with one additional 
parameter through the weight function 

T(p + l)z x 

u{x) = T(p + x + i)T( x + iy z > °' 

with (3 a parameter such that (3 > — 1. This gives the discrete weight 

w( *> = crnpi' z > °- (4 - 3) 

where (/3 + 1)& = T(/3 + 1 + k)/T(/3 + 1) is the Pochhammer symbol, on the lattice N. 
The weight (4.3) satisfies the discrete Pearson equation (3.10) with 

a(k) = k(k + p), r{k) = -k 2 - f3k + z, 

and so the generalized Charlier polynomials are semi-classical orthogonal polynomials. 
The special case /3 — was first considered by Hounkonnou, Hounga and Ronveaux [33] 
and later studied by van Assche and Foupouagnigni [61]. 

For the generalized Charlier weight (4.3), the orthonormal polynomials p n (k; z) 
satisfy the orthogonality condition 



^^Pmik; z)p n (k; z )jw~y: r , — ^m,n, 



k=0 v ' 

and the three-term recurrence relation 

xp n (x; z) = a n+1 (z)p n+1 (x; z) + b n (z)p n (x; z) + a n (z)p n ^ 1 (x; z), (4.4) 

with pi(x; z) — and po(x; z) = 1. Our interest is determining explicit expressions for 
the coefficients a n (z) and b n (z) in the recurrence relation (4.4). 

Smet and van Assche [56, Theorem 2.1] proved the following theorem for recurrence 
coefficients associated with the generalized Charlier weight (4.3). 

Theorem 4.1. The recurrence coefficients a n (z) andb n (z) for orthonormal polynomials 
associated with the generalized Charlier weight (4-3) on the lattice N satisfy the discrete 
system 

( a l+i - z )( a l ~z) = z(b n - n)(b n -n + (3), 

(4.5) 

b n + - n + {3 + 1 = nzj a 2 n1 
with initial conditions 

4 = 0, b = ^^f ] = ^ ln ^ m W^)X (4-6) 
with I u (x) the modified Bessel function. 
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Smet and van Assche [56, Theorem 2.1] show that the system (4.5) is a limiting 
case of a discrete Painleve equation, namely the first dPpy in [60, p. 723]. Since the 
initial conditions (4.6) involve modified Bessel functions, then clearly the solutions of 
the discrete system (4.5) are expressed in terms modified Bessel functions. 

Using the discrete system (4.5) and the Toda system (3.8), Filipuk and van Assche 
[20] show that the recurrence coefficient b n can be expressed in terms of solutions of a 
special case of P v which can be transformed into Pm- However their proof is rather 
involved and several details are omitted due to the size of expressions involved. Further 
Filipuk and van Assche [20] do not obtain explicit expressions for the the recurrence 
coefficients a n and b n . 

The relationship between the recurrence coefficients a n (z) and b n (z) and classical 
solutions of Pur can be shown in much more straightforward way and also obtain explicit 
expressions for these coefficients. First we obtain explicit expressions for the moment 
/jio(z) and the Hankel determinant A n (z). 

Theorem 4.2. For the generalized Charlier weight (4-3) the moment fJ>o(z) is given by 

OO f. 

fc=0 ^ > k 

with I v (x) the modified Bessel function, and the Hankel determinant A n (z; (3) is given 
by 

A n (z;P) = [r(/3 + l)]X(^ /2 /,(2^)). (4.8) 
Proof. Since the modified Bessel function I u (x) has the series expansion [53, §10.25.2] 

oo 



o k\r(u + k + iy 

then the expression (4.7) for the moment (J,o(z) follows immediately. Then using 
Theorem 3.1, we obtain the expression (4.8) for the Hankel determinant A n (z; (5). □ 

Hence we obtain explicit expressions for the recurrence coefficients a n (z) and b n (z). 

Corollary 4.3. The coefficients a n (z) and b n (z) in the recurrence relation (4-4) have 
the form 

<4« = (4) 2 (,„a„(^)), „, M - 4 (>»^f ) • («) 

with A n (z; (3) given by (4-8). 

Proof. These follow immediately from Theorem 3.3. □ 

Finally we relate the Hankel determinant A n (z; /3) to solutions of (2.12), the Pm' 
cr-equation. 
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Theorem 4.4. The function 

S n (z;f3) = z-^\nA n (z;/3), 
with A n (z; (3) given by (4-8), satisfies the second- oder, second-degree equation 



(4.10) 



d 2 S n 
' dz 2 



n — (n + (3) 
. dS„ ( dS, 



dSn 

dz 



- 1 



1 2 



z—^- - S n + \n(n - 1) 



(4.11) 



dz \ dz 

Proof. Equation (4.11) is equivalent to (2.12) through the transformation 

S n (z; P) = a(z) + \z+ \n 2 - \ n {fi + 1) - ±/3 2 , (4.12) 

with parameters (-^o, ^oo) = {n + (3, n — (3), as is easily verified. Then comparing (4.12), 
with S n given by (4.10), to (2.15), with v — f3, e x — — 1 and e 2 = 1, gives the result. □ 

Remarks 4.5. 

(i) In terms of S n (z; (3) given by (4.10), the coefficients a n (z) and b n (z) in the recurrence 
relation (4.4) have the form 

a 2 n (z) = z^S n (z; 0), b n (z) = S n+l (z; 0) - S n (z; fi). 

(ii) Solving (4.12) for a(z) and substituting in (2.13) yields the solution of the 
Hamiltonian system (2.9) given by 

zS' 1 i(z)-2(n + /3)S' n (z)+2n 



q(z) ±S' n (z)[l - S' n (z)} 

with parameters (i?o, ^oo) = (n + (3,n — (3). 



5. Meixner polynomials and generalizations 

5.1. Meixner polynomials 

The Meixner polynomials M n (k;a,z) are a family of discrete orthogonal polynomials 
introduced in 1934 by Meixner [44] given by 

M n (k; a, z) — 2 F 1 (-n, —k; -a; 1 - 1/z) , < z < 1, (5.1) 

with a > 0, where 2 F\ (a, b; c; z) is the hypergeometric function, see [4, 14, 35, 53]. In 
the case when a = —N, with iV e N and z — p/(l — p), these polynomials are referred 
to as the Krawtchouk polynomials for k e {0, 1, ... , iV} 

K(k;p,N) = 2 F 1 (-n,-k;N;l/p), < p < 1, (5.2) 

which were introduced in 1929 by Krawtchouk [39]. The Meixner polynomials (5.1) are 
orthogonal with respect to the discrete weight 

w (jfe) = ^ki! > tt >0, z>0, (5.3) 
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and satisfy the orthogonality condition 

°° (a)kZ k n\ z~ n 

M m{k] a, z)M n (k; a, z) — - — = Wl _ x Q <W 

k=0 K - l a JnU Z) 

The weight (5.3) satisfies the discrete Pearson equation (3.10) with 

a(k) = k, r(k) — (z — l)k + za. 
From (5.3), the moment fJ>o( z ) is given by 



00 



k=0 

Hence from Theorem 3.1 that the Hankel determinant A n (z) is given by 

n(n-l)/2 

A n (z) = wm = _ (w+a _ 1} U k ^+ k r~ k -\ 

^ > k=l 

and so, from Theorem 3.3, the recurrence coefficients are given by 
9 . . n(n + a — l)z , , , n+(n + a)z 

4W = (1 _, )2 ■ M«> = ■ 

5.2. Generalized Meixner polynomials 

In a similar way to that for the Charlier weight above, Smet and van Assche [56] 
generalized the Meixner weight (4.2) with one additional parameter through the weight 
function the weight function 

. T{a + x)T{(5)z x 

u{x) = r(n)r ( ;.,)n,. ,r z > °« 

with a, j3 > 0, which gives the weight 

wW = K 2 > a (5 ' 4) 

The weight (5.4) satisfies the discrete Pearson equation (3.10) with 

a(k) = k(k + (3- 1), r(k) = -k 2 + (z + 1 - (5)k + za, 

and so the generalized Meixner polynomials are semi-classical orthogonal polynomials. 

Boelen, Filipuk and van Assche [5] considered the special case of (5.4) when (3 = 1 
and showed that the recurrence coefficients a n and b n satisfy a limiting case of an 
asymmetric dPpv equation. We note that the special case a = (5 gives the classical 
Charlier weight (4.2) and the case a = 1 corresponds to the classical Charlier weight on 
the lattice N + 1 - /3. 

For the generalized Meixner weight (5.4), the orthonormal polynomials p n (k; z) 
satisfy the orthogonality condition 



k 
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and the three-term recurrence relation 

xp n (x; z) = a n+1 (z)p n+l (x; z) + b n (z)p n (x; z) + a n (z)p n _ 1 (x; z), (5.5) 

with pi(x;z) = and po(x; z) = 1. As for the generalized Charlier weight (4.3), our 
interest is determining explicit expressions for the coefficients a n (z) and b n (z) in the 
recurrence relation (5.5). 

Smet and van Assche [56, Theorem 2.1] proved the following theorem for recurrence 
coefficients associated with the generalized Meixner weight (5.4). 

Theorem 5.1. The recurrence coefficients a n (z) andb n {z) for orthonormal polynomials 
associated with the generalized Meixner weight (5.4) on the lattice N satisfy 

a 2 n = nz — (a — l)x n , 

b n = n + a — (3 + z — (a — l)y n /z. 

where x n and y n satisfy the discrete system 



i \i , a — 1 . ( a — (3\ 

[Xn + yn){Xn+l + Vn) = — 2/n(t/n ~ Z) \ V n ~ Z 1 

(a - l)x n {x n + z) ( a- /3 

[Xn + y n ){Xn + Vn-l) = —, T~\ [X n + Z 



(5.6) 

(a — l)x n — nz \ n a — 1 / 
with initial conditions 

i n 7 az M(a + 1, j3 + 1, z) d n „ w n . . 

and M(a,j3,z) is the Kummer function. 

We note that M(a,/3,z) = \Fi(a, f3, z), the confluent hypergeometric function [53, 
§13]. Smet and van Assche [56] show that the discrete system (5.6) can be identified as 
a limiting case of an asymmetric dPrv equation. Filipuk and van Assche [19] show that 
the system (5.6) can be obtained from a Backlund transformation of Py (1-2). Since 
the initial conditions (5.7) involve Kummer functions, then clearly the solutions of the 
discrete system (5.6) are expressed in terms Kummer functions. 

Using the discrete system (5.6) and the Toda system (3.8), Filipuk and van Assche 
[19] show that the recurrence coefficients a n (z) and b n (z) are related to classical solutions 
of P v (1-2), for the parameters 

(A, B, C, D) = {\{a- l) 2 , -\{n + a - /3) 2 , ±{n + /3), -\). 

However their proof is rather cumbersome and most of the details are omitted due to the 
size of expressions involved. Further Filipuk and van Assche [19] do not obtain explicit 
expressions for the the recurrence coefficients a n and b n . 

In an analogous way to that for the generalized Charlier polynomials in §4.2 
above, the relationship between the recurrence coefficients a n (z) and b n (z) and classical 
solutions of Py can be shown using a much more straightforward way and also obtain 
explicit expressions for these coefficients. First we obtain explicit expressions for the 
moment [J,o(z) and the Hankel determinant A n (z; a, (3). 
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Theorem 5.2. For the generalized Meixner weight (5.4) the moment Ho{z) is given by 

(5.8) 



£fefe = M(a,/M, 



k=0 



(f3) k k\ 



with M(a,(3,z) the Kummer function, and the Hankel determinant A n {z) given by 

A n (z;a,f3)=W n (M(a,f3 : z)). (5.9) 

Proof. Since the Kummer function M(a, f3, z) has the series expansion [53, §13.2.2] 



M(a,(3,z) = J2 



a)k k 

z i 



k=0 



(P)kk\ 



then the expression (5.8) for the moment (J®(z) follows immediately. Then we use 
Theorem 3.1, to obtain the expression (5.9) for the Hankel determinant A n (z; a, (3). □ 

Hence we obtain explicit expressions for the recurrence coefficients a n (z) and b n (z). 

Theorem 5.3. The coefficients a n (z) and b n (z) in the recurrence relation (5.5) have 
the form 



a n(*) = (mA n (z;a,/3)), b n (z) = z- 

with A n (z) given by (5.9). 



In 



z;a,(3) 



dz \ A n (z;a,P) 



Proof. These follow immediately from Theorem 3.3. 



(5.10) 



□ 



Finally we relate the Hankel determinant A n (z) to solutions of (2.23), the P v o- 
equation. 

Theorem 5.4. The function 



S n (z; a, (3) = z— In A n (z; a, /3), 
dz 

with A n (z;a,{3) given by (5.9), satisfies the second-oder, second-degree equation 



(5.11) 



d 2 S n 
' dz 2 



[z + n + (5- 1)— - -S n - \n{n - 1 + 2a) 
dz 1 



-4- 



, dS„ ( dS, 



dz \ dz 



— n — a + (3 



dSn 

' dz 



-S n + \n{n - 1) 



(5.12) 



Proof. Equation (5.12) is equivalent to (2.23) through the transformation 

S n (z; a, /3) = a(z) + \{2a - (3 + 3n - l)z 

+ |n 2 + \{2a -3(3- l)n + §(2a - (3 - if 



(5.13) 



with parameters 

«i = i(2a-/3-n-l), « 3 = \{2a - j3 + 3n - 1), 

k 2 = -\{2a + /3 + n - 3), k 4 = -\{2a - 3/3 + n + 1), 

as is easily verified. Then comparing (5.13), with S n given by (5.11), to (2.27a) gives 



(5.14) 



the result. 



□ 
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Remarks 5.5. 

(i) In terms of S n (z;a,j3) given by (5.11), the coefficients a n (z) and b n (z) in the 
recurrence relation (5.5) have the form 



(hi) Solving (5.13) for a(z) and substituting in (2.25) yields the solution of the 
Hamiltonian system (2.21) given by 



for the parameters (5.15). 
6. Discussion 

In this paper we have studied semi-classical generalizations of the Charlier polynomials 
and the Meixner polynomials. These discrete orthogonal polynomials satisfy three- 
term recurrence relations whose coefficients depend on a parameter. We have shown 
that the coefficients in these recurrence relations can be explicitly expressed in terms 
of Wronskians of modified Bessel functions and Kummer functions, respectively. 
These Wronskians also arise in the description of special function solutions of the 
third and fifth Painleve equations and the second-order, second-degree equations 
satisfied by the associated Hamiltonian functions. The results in this paper are more 
comprehensive than those in [20] for generalized Charlier polynomials and in [5, 19] 
for generalized Meixner polynomials. The link between the semi-classical discrete 
orthogonal polynomials and the special function solutions of the Painleve equations 
is the moment for the associated weight which enables the Hankel determinant to be 
written as a Wronskian. In our opinion, this illustrates the increasing significance of the 
Painleve equations in the field of orthogonal polynomials and special functions. 
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a 2 n (z) = z—S n (z; a, fj), b n (z) = S n+1 (z; a, (3) - S n (z; a, (3). 
(ii) Substituting the parameters (5.14) in (2.24) yields 
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(5.15) 



p(z) 




2zS'^ + 4[S' n ] 2 - 2{z-2a + P-3n + l)S' n + 2S n + n{n + 2a - 1) 

4[S' n - a - n + l][S' n - n] 
2zSl - 4[S' n } 2 + 2 (z - 2a + g - 3n + 1) S' n - 2S n - n(n + 2a - 1) 

4[S' n -a + f3- n] 
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